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ABSTRACT 

The recently discovered J0651+2844 is a detached, eclipsing white dwarf binary with an orbital 
period of 765 s. We investigate the prospects for the detection of gravitational radiation from this 
system and estimate the effect of the tidal deformation of the low- mass component on the period evo- 
lution of the system. Because of the high inclination of the system, the amplitude of the gravitational 
waves at Earth will be as much as a factor of two lower than that from an optimally oriented sys- 
tem. The dominant contribution of tidal corrections to the period evolution comes from the increase 
in rotational energy of the components as they spin up to remain tied to the orbital period. This 
contribution results in an advance of the timing of the eclipses by an additional 0.3 s after one year. 
Subject headings: Binaries: eclipsing — Gravitational waves — White dwarfs 
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1. INTRODUCTION 

The Galactic population of white dwarf bi- 
naries is expected to produce a foreground 
of gravitational radiation signals in the milli- 
hertz frequency band fevans. Iben fc Smarr 1 19871 : 
Hils, Bender & Webbink' 1990; Hils & Bender '200C| 
Nelemans 2003; Benacquista, DcGocs & Lunder 2004 
Benacquis ta fc HoUev-Bockelmannl 120061 : IRuiter et ahl 
2010). Most of the known binaries that will be de- 
tectable by planned space-based detectors such as 
the Laser Interferometer Space Antenna (LISA) are 
interacting systems, and so the period evolution of these 
systems will be dominated by mass transfer. The recent 
discovery of SDSS J065133.33-h284423.3 (hereafter 
J0651), a detached eclipsin g white dwarf bina ry with a 
12 minute orbital period (jBrown et al.l 120111 ) provides 
the opportunity for indirect observation of gravitational 
radiation through the measurement of the period evolu- 
tion. If detectable by a future detector such as LISA, 
the combination of direct measurement of gravitational 
radiation and the observation of the period evolution 
through timing of the eclipses in the intervening years 
will allow for detailed study of the internal structure 
of the component white dwarfs through the dissipation 
of orbital energy via tidal deformation. In this paper, 
we evaluate the expected contribution of the tidal and 
rotational deformation of the components to the period 
evolution of the binary, and the effect on the direct 
measurement of gravitational radiation with future 
detectors. 

J0651 is an eclipsing binary with an angle of inclina- 
tion of i = 86.9ti-^ degrees (jBrown et al.l 1201 If ). They 
report an orbital period of P = 765.2062 ± 0.003 s 
for J0651, with component masses of Mp — 0.25 Mq 
and AIs — 0.55 Mq. Ellipsoidal variations in the light 
curve indicate a 5% difference in the projected area 
of the primary, showing that it is tidally distorted. 
The sinusoidal velocity curve indicates that the orbit 
is circular, while the symmetric appearance of the 
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ellipsoidal variations indicates that the primary is nearly 
tidally locked. Both the primary and secondary have 
radii that are slightly larger than those predicted by 
models of isolated white dwarfs. The primary has an 
observed radius of 0.0353 ± 0.0004 Rg, which is 5% 
larger than models predict (|Panei et all [20071) . The 
secondary has an observed radius of 0.0132 ± 0.0003 Rq, 
which is 4% greater than models predict. These 
variations in the sizes and shapes of the components 
indicate that they are tidally distorted due to their close 
proximity. This distortion will have an effect on the 
orbital period evolution. If this effect is not properly 
taken into account, calculated values of the binary 
properties from gravitational radiation observations 
will be in error. Previous studies of the effect of tidal 
energy dissipation and apsidal precession on detached 
binaries have focused on eccentric syste ms, where the 
apsid al precession can be measured (IValsecchi et all 
[20T1I: IWillems. Vecchio fc Kaloeeral I2008D or 
where the va rying tidal stresses excit e inter- 
nal modes (j'Willems. Delove fc Kalogeral iOlO), 
or circul ar systems w ith non-rotating compo- 
nents (.Fuller fc Laill2011afl. or through ele ctromagnetic 
interactions (jWu. Willes fc: Ramsavl I2008D . AUhough 
the components of J0651 appear to be tidally locked, 
the gradual decay of the orbit due to gravitational 
radiation requires there to be a slight mismatch between 
the spin frequency of the components and the orbital 
frequency of the binary. The slow spin up of the 
components as a result of the torques applied due to this 
mismatch will provide the bulk of the deviations from 
pure gravitational wave emission in the period evolution 
of the binary. 



2. ORBITAL EVOLUTION 

The tidal distortions produce non-zero quadrupole mo- 
ments in the components. If we include the quadrupole 
moments in the expansion of the potential, we find: 
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where M is the total mass of the system, ^ = MpMg/M 
is the reduced mass, and Q is the quadrupole moment. 
Here, we have assumed that the secondary is not signif- 
icantly distorted (as will be shown later), and that the 
quadrupole moment of the primary is symmetric so that 
it can be described by a single value. This potential re- 
sults in a variation of Kepler's law, so that the orbital 
angular frequency (w) and orbital radius (a) are related 
by: 

Since the orbital radius is not directly measurable in 
gravitational radiation, the orbital angular frequency is 
the variable of interest. We will show later that the 
quadrupole correction will be less than 1%, so we can 
solve for a as a perturbation about the Keplerian solu- 
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tion, oo = (GM/w^) . We set a = ao{l + Sa/ao), and 
find the perturbation to be: 
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We are now in a position to calculate the gravitational 
wave strain and the gravitational wave power for a per- 
turbed system. Assuming point masses for the orbital 
configuration in order to calculate the total quadrupole 
m oment of the binary sys tem following the method 
of IPeters fc Mathewd (|l96l . we find that the power ra- 
diated through gravitational radiation is: 
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The energy radiated by gravitational waves comes from 
the orbital energy of the system, and this manifests it- 
self in a gradual shrinking of the orbit with an increase 
in the orbital frequency, known as the "chirp." In the 
case of J0651, the tidal distortion of the components pro- 
vides a torque that spins up the components so that their 
rotations remain locked to the orbital period. The in- 
crease in the rotational kinetic energy of the components 
also comes from the orbital energy of the system, and so 
the increase in orbital frequency will be larger than that 
solely due to the emission of gravitational radiation. 

If the torques on the tidal deformations couple to any 
internal modes for each component, additional energy 
can be lost to heating of the component stars. However, 
the resonance condition for coupling the tidal deforma- 
tions to these internal modes requires that the driving 
frequency of the distortion be comparable to the natural 
frequency of one of the modes. The driving frequency is 
given by: 

cjf 2 (w - r?) (5) 

where uj is the orbital freq uency and fl is the s pin fre- 
quency of the white dwarf (jFuUer fc Laill2011 bV bmce 
the lowest normal mode frequency in a whit e dwarf is on 
the order of a mHz (jWinget fc Keplerll2008[ ). and the or- 
bital frequency is 1.3 mHz, then < uj/2 in order for ujt 
to be comparable to any normal mode frequency. This 
would be clearly visible in the light curve, and so we can 
neglect coupling to the internal modes. An additional 



channel for energy loss in the system may be through 
heating without any direct coupling to normal modes. 
[Fuller fc Lail (|2011bl ) have shown that the tidal heating 
rate can be compared to the tidal spin-up energy rate by: 
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and so if ~ uj, then the energy loss due to heating 
can also be neglected. A more detailed treatment of the 
effects of tidal heating on the spin evolution of J0651 
and the implications of a non-negli gible contribution to 
heating can be found in lPirol (pOllI ). 

Here, we obtain a lower bound on the increase in the 
chirp by only considering the rotational kinetic energy 
and ignoring any internal heating due to coupling with 
internal modes, thus the change in total energy of the 
system is equal to the power radiated through gravita- 
tional radiation, so: 
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where Ip and Is are the moments of inertia of the primary 
and the secondary, respectively. Introducing the dimen- 
sionless quantities due to the quadrupole correction to 
the potential (Aq) and the spin- up of the components 
(Ai), where: 
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where ujq is the value of the chirp in the absence of any 
tidal deformations: 



we find 
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The chirp can be used in combination with the grav- 
itational wave frequency (wqw = 2w) to determine the 
"chirp mass" , a combination of the reduced mass and the 
total mass given by: ^Ac — p^^^M"^!^ . If the corrections 
due to tidal deformation are not properly accounted for, 
the calculated chirp mass (TWc) will be greater than the 
true chirp mass by: 
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The amplitude of the two polarizations are given 
by fe.g.: ICutleii[T998h as: 
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where d is the distance to the binary and i is the in- 
clination. We note that the fact that J0651 is eclipsing 
implies cosi ~ 0.05 and so the x -polarization ampli- 
tude is near zero and the -|— polarization amplitude is 
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near its minimum. Thus, the strain amphtude may be 
as much as a factor of two lower than the value estimated 
in iBrown et alj ()201lD . Taking into account the pertur- 
bation in a due to the deformation of the components 
(see Equation Is]) , we can express the basic amplitude as: 
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Using the chirp mass calculated from oj and Cj, to elim- 
inate ^J^M'^/^ from /i, iSchutzl ()1997[) showed that the dis- 
tance to the binary can be calculated using: 
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However, this expression is only valid for point masses 
and assumes that the chirp is entirely due to energy loss 
through gravitational radiation. Including effects due 
to the tidal and rotational distortions of the component 
masses results in the following equation for the distance: 

Thus, we see that including effects due to tidal defor- 
mation yields corrections to Kepler's law, and calculated 
values of the chirp, chirp mass, and distance to the bi- 
nary. All of these corrections can be expressed in terms 
of Aq and Aj. In the next section, we estimate the size 
of these corrections in J0651. 

3. PERTURBATIONS IN J0651 

Both components of J0651 are relatively cool white 
dwarfs with masses substantially below the Chan- 
drasekhar mass. Conseq uently, we can conside r them 
to be n = 1.5 polytropes. ISirotkin fc KimI ()2009l ) looked 
at perturbations in the internal structure for polytropes 
in close binaries using a self-consistent field method to 
calculate the effects of rotation and tidal deformation in 
synchronized, circular-orbit binaries. Their calculations 
assume that one component of the binary is a point mass 
and they then determine the perturbations of the other 
component. The volume radius R is defined as the radius 
of a sphere with the equivalent volume of the distorted 
polytrope, and the dimensionless volume radius (^i) is 
determined by: 



where 



^1/3 



1/2 



(18) 



(19) 



and K and pc are the pressure constant and the central 
density, respectively. The mean density is determined by 
the mass (M*) of the star and the volume radius, so that 
p — (47ri?'^). A dimensionless description of the 

mass concentration is given by Pc/ P- The moment of in- 
ertia (/), is characterized by the dimensionless quantity: 

1 



a^Pc 



p{a 



(20) 



where the orbital plan e is assumed to lie in the xy-plane. 
ISirotkin fc KimI ()2009[ ) produce numerical fitting formu- 
lae for ^1, /, and Pc/p in terms of the relative radius 



TABLE 1 

Structure parameters of ,10651 
components and unperturbed n = 1.5 
polytropes. 





Pc/p 




/ 


Primary 


6.0486 


3.6783 


92.0131 


Secondary 


5.9901 


3.6538 


90.9334 


Unperturbed 


5.9907 


3.6538 


93.156 



Ra ~ R/a. We can recover the moment of inertia from 
the dimensionless quantities, using Ra, a, and M^: 



(21) 



Since any corrections to the orbital separations will result 
in second order corrections to the orbital evolution, we 
assume that the orbital separation is the unperturbed 
separation for these calculations, and so a — 0.1676 Rq. 

For the primary, we have Mp = 0.25 Mq, Ra = 
0.21, and the mass fr action M„/M = 0.3. Using Ta- 
bles 4, 5, and 7 from ISirotkin fc KimI H2009f ). we find 
Pc/p = 6.0486, = 3.6783, and / = 92.0131. This 
corresponds to a moment of inertia for the primary 
of Ip = 6.111 X 10"^ MqRI = 5.912 X lO''^ kg-m^. 
For the secondary, we have Ms = 0.55 M©, Ra = 
0.08, and the mass fraction Mg/M — 0.7. This gives 
Pc/p = 5.9901, 6 = 3.6538, and / = 90.9334, corre- 
sponding to a moment of inertia for the secondary of 
/, = 1.914 X 10"^ MqR2, = 1.851 X lO'*^ kg • m^. These 
values, along with those for an unperturbed, isolated 
n = 1.5 polytrope are given in Table [TJ From here we 
can see that the secondary is nearly identical to an un- 
perturbed white dwarf with the exception of a lower mo- 
ment of inertia. We can interpret this as implying that 
the secondary is relatively undistorted. Thus, we need 
only consider the quadrupole moment of the primary. 

In order to estimate the quadrupole moment of the 
primary, we note that the quadrupole moment is pro- 
portional to the moment of inertia with a proportion- 
ality constant related to the structure of the star. In 
this case, we do not have a detailed description of the 
structure of the primary, but we can estimate an upper 
bound on the quadrupole moment by assuming constant 
density and a shape that deviates from spherical by a 
small amount. Since the light curve of J0651 indicates 
a 5% ellipsoidal variation of the primary. Thus, we as- 
sume a 5% deviation from spherical and compute the 
quadrupole moment. It can be shown that to first or- 
der in the deviation, the quadrupole moment is 5% of 
the moment of inertia. Therefore, we take as an upper 
bound, Q = 0.05/p = 2.956 x 10^^ kg . ni2. 

Taking the above values for /p, /s, and Q and using the 
observed values of the masses and frequency, we find that 
the dimensionless perturbations of Equations [5] and [S] 
are: 



Aq = 1.36 X 10"^ 
Ai = 0.0166. 

4. DISCUSSION 



(22) 
(23) 



The perturbations introduced by the distortion and 
tidal locking of the components of J0651 will be on the 
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order of a few percent. Nearly all of this will come from 
the spin- up of the components due to the torques induced 
through the inspiral in order to maintain tidal locking. 
The effects of a non-zero quadrupole moment will con- 
tribute less than 0.1% to any perturbation, and can safely 
be ignored at the separations found in J0651. If the effect 
of spin-up is ignored in future gravitational observations, 
then the chirp mass will be over-estimated by 3%, which 
is well within the errors expected from spectrum mod- 
eling. The distances will then also be over-estimated if 
the standard formula of Schutz (19911 is used. The cal- 
culated distance will be nearly 5% greater than the true 
distance. At a current estimated dis tance of 1 kpc, an d 
an apparent magnitude of g = 19.1 ([Brown "eraLl[20TTI) . 
GAIA can be expected to provide a distance accuracy 
of around 8-20% ( Bail er- Jonei 120091 ) . Thus, ignoring 
the perturbations due to tidal distortions and tidal lock- 
ing when using gravitational wave observations will still 
provide better estimates of the distance and masses of 
white dwarf binaries such as J0651 than can be obtained 
through electromagnetic observations. 

One area in which electromagnetic observations will 
be superior to gravitational wave observations will be 
in the measurement of the chirp. An estimate of the 
minimum measurable chirp through gravitational wave 
observations can be found by requiring the gravitational 
wave frequency of the binary to shift by one resolvable 
frequency bin during the observation time Tobs- Since the 
size of the resolvable frequency bins scales as the inverse 
of the observation time, the minimum measurable chirp 
is estimated by /gw ^ ^obs' where /gw — uj/tt. Ignoring 
the perturbations, w = 9.184 x 10^^^ Hz^, which implies 
a 6 year observation time. The large signal-to-noise ratio 
expected from this source, combined with improved data 
analysis may shorten this requirement to two or three 
years. On the other hand, timing of the eclipses allows 
for the chirp to measured within one year. After a time 
interval T, the eclipses will occur earlier by an amount: 



If the perturbations are ignored, St — 5.57 s for T = 1 yr. 
Including the perturbations shifts this to dt = 5.85 s — 
a difference of less than 0.3 s. However, after two years, 
this difference grows to more than 1 s. 

Combining both electromagnetic and gravitational 
wave observations will allow for the properties of J0651 
to be measured to better precision than either observa- 
tion alone can achieve. Assuming that a space-based 
observatory such as LISA is launched in the early 2020's, 
the timing of the eclipses of J0651 will have measured 
an accumulated advance of the eclipse of nearly 10 min- 
utes, 30 seconds of which can be attributed to tidal ef- 
fects. Thus, estimates of the distance to J0651 can be 
made using the gravitational wave chirp as measured by 
eclipses and the strain amplitude as measured by LISA 
without having to wait for the chirp to be measured in 
the gravitational wave signal. Furthermore, astrometric 
measurements of J0651s distance will have been made by 



GAIA, and this can be compared with the gravitational 
wave distance. Keeping in mind that the perturbations 
due to the tidal interactions that we have calculated here 
ignore any thermal dissipation of energy into heating of 
the white dwarfs, the true perturbation may be large 
enough to cause a discrepancy between the gravitational 
wave distance and the astrometric distance. 

5. CONCLUSION 

We have shown that the tidal deformation of the com- 
ponents of J0651 will cause the orbital evolution to devi- 
ate from that predicted by the simple loss of energy due 
to the emission of gravitational radiation as some orbital 
energy will go into spinning up the components in or- 
der to maintain tidal locking. Ignoring any additional 
dissipation of energy through heating of the component 
white dwarfs, this extra energy sink will cause the chirp 
to be about 5% greater than due to gravitational radia- 
tion alone. This is a lower bound on the perturbation, as 
thermal dissipation will only extract more energy from 
the orbit. Although the deformation will also cause the 
orbit to deviate from pure Keplerian, we have shown that 
this effect will contribute less than 0.1% to the chirp, and 
can probably be safely ignored. 

As J0651 is an eclipsing system, its inclination is quite 
high. This has a disadvantage for gravitational wave ob- 
servations, since the gravitational wave strain is more 
than a factor of two weaker for this orientation com- 
pared with the optimal orientation of face-on. On the 
other hand, timing of the eclipses allows for the chirp 
to measured to better accuracy over a shorter observing 
time than can be done with gravitational wave observa- 
tions. Since the precision of this measurement increases 
with observing time and no gravitational wave observa- 
tory in this frequency band is expected for at least a 
decade, eclipse timing will provide an exquisitely mea- 
sured chirp for use with gravitational wave observations 
as soon as the detectors come online. In the intervening 
years between now and the launch of millihertz gravita- 
tional wave observatories, GAIA will have completed its 
mission and analyzed parallax measurements for J0651 
(if it is targeted). The distance obtained through paral- 
lax is expected to have errors of 8-20%, while the great- 
est errors in the distance obtained through gravitational 
waves will be due to the uncertainties in the amount of 
energy dissipated through tidal interactions. We have 
estimated this to be smaller than the errors in parallax 
measurements. Thus, millihertz gravitational wave ob- 
servations will provide better measurements of the dis- 
tances to detached double white dwarfs such as J0651 
and other faint objects than can be obtained through 
parallax measurements, with the dominant source of er- 
ror being the tidal interactions. 
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